Gravitational Lamb Shift of Bose-Einstein Condensates due to Spacetime Fluctuations 
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We investigate the oscillation of the center of mass of trapped Bose-Einstein condensates coupled 
to the zero-point fluctuations of the gravitational field. A semiclassical analysis is performed that 
allows to calculate the mean square amplitude of the oscillation. In analogy with the Lamb shift 
in quantum electrodynamics, this gives rise to an upshift of the energy of the trapped condensates. 
We show that for an elongated trap, the energy shift scales quadratically with the length as well 
as cubically with the total number of atoms, leading to an energy increase of 1 % using a 5 cm 
long trap with 10* rubidium atoms. This could potentially lead to the first observable effect of low 
energy quantum gravity and provide a stringent test for any viable theory of quantum gravity. 
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Introduction and summary. — Physics at the turn of the 
21st century has seen tremendous advances in confirming 
and controlling the quantum behavior of complex mat- 
ter [ll-Q . The unprecedented precision of atom interfer- 
ometry holds out new prospects for laboratory tests of 
general relativity [1] . The development of large molecule 
interferometry is providing an arena to challenge our un- 
derstanding of the foundations of quantum mechanics Q . 
Furthermore, there has been a surge of interest in detect- 
ing signatures of low energy quantum gravity and unified 
theories using quantum matter waves [6|-|8|. 

The lowest energy quantum gravity effect one would 
hope to probe is the zero-point vacuum fluctuations of 
spacetime. A possible scenario is the decoherence of mat- 
ter waves through coupling to a fluctuating metric 
Unfortunately, such efforts have been hampered by the 
smallness of the effect and hence high sensitivity require- 
ment [13], the need to suppress other environmental de- 
coherence effects within the interferometer Q, and the 
ambiguity of predictions in the absence of a consistent 
quantum theory of gravity [l3 |. 

Skeptics of quantum gravity phenomenology often ar- 
gue that its experiments require accessing energies up to 
the Planck scale of 10^^ GeV, which is far unattainable 
from our low energy environment. However, this view 
must be assessed with caution. After all, quantum elec- 
trodynamics (QED) as a highly successful quantum field 
theory, has demonstrated how renormalization involv- 
ing high energy states can influence low energy physical 
observables. The Casimir effect testifies to the observ- 
able nature of the modified renormalization of quantum 
vacuum through boundary conditions set by conducting 
plates at low energy. 

It is interesting to note that the decoherence of elec- 
trons due to vacuum fluctuations of the electro mag netic 



(EM) field has been predicted for some time [15|, but 
is still far below measurable levels One important 

consequence of the vacuum fluctuations of a quantum 



field is the energy shift of a bound state coupled to the 
fleld. The measurement of the energy or corresponding 
frequency shift through spectroscopy or resonance tech- 
niques can generally be done more accurately than the 
visibility measurement for quantum decoherence. Indeed, 
the resulting Lamb shift |17| of electron energy levels in 
an atom provides the first experimental evidence for EM 
vacuum fluctuations and agrees with predictions from the 
renormalization procedures in QED. 

In this paper we investigate the energy shift of quan- 
tum bound states due to spacetime fluctuations. We 
shall refer to this phenomenon as the "gravitational Lamb 
shift" , even though the quantum states may be conflned 
by non-gravitational interactions. The possibility of en- 
ergy shift due to spacetime fluctuations has been pre- 
viously noticed [18]. Here we proposed a realistic sce- 
nario where such an effect could be measured. We have 
in mind a large number of trapped cold atom and but 
trapped molecules may also be relevant. For such a sys- 
tem, the gravitational Lamb shift is most important for 
its center-of-mass (CM) wavefunction. However the size 
of the effect is generally quite small. 

A natural analogy with gravitational wave detectors 
fortuitously emerges in this work, which makes possible 
significant amplifications of the effect by elongated coher- 
ent atom traps. The trap acts as a quantum detector for 
zero-point energy gravitons. The resulting gravitational 
Lamb shift not only scales cubically with the number of 
atoms, but also quadratically with the longitudinal ex- 
tension of the trap, leading to an energy increase of 1 % 
using a 5 cm long trap with 10® rubidium atoms. The 
required number of cold atoms and the physical size of 
the atom trap are not far-fetched from currently being 
achieved in the laboratory 1^ 2o| . 

Position oscillations under vacuum EM fluctuations. — 
Let us begin by recapitulating Welton's theory of the 
Lamb shift of a non-relativist electron in a hydrogen atom 
2lj . The phenomenon has a neat semiclassical explana- 
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FIG. 1: In the synchronous coordinates {t,x,y,z), trajecto- 
ries with constant x,y,z are freely falling. For hij — 0, these 
represent stationary particle positions. With hij 7^ 0, a par- 
ticular wave component, say in the z-direction, modulates the 
distance between any point on the x-y plane and a chosen ori- 
gin along a macroscopic straight line, defined as those on the 
averaged flat space. For a particular polarization the instan- 
taneous displacements on macroscopic concentric circles are 
illustrated in the diagram on the left. As in the standard 
description of weak GWs, the radial displacements are pro- 
portional to the radius. The position oscillations along e.g. 
the a;-axis is illustrated in the diagram on the right. Eq. (j6]) 
is obtained by integrating the resulting mean square devia- 
tions in all wave directions. In both diagrams the dashed 
lines represent displaced positions relative to the origin. 



tion in terms of a small gain of the electron's average 
potential energy due to its position oscillation induced 
by a stochastic classical EM radiation field at zero tem- 
perature T = 0. 

The EM fiuctuations at a finite T are homogeneous 
and isotropic vifith a mean square electric field given by 
the statistical average 
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(1) 



in terms of the thermal energy distribution of a harmonic 
oscillator with angular frequency uj: 



E{lu,T) = 



huj 



hid 



2 exp{hLj/kT)- 1' 



(2) 



Relations ([TJ and 1^ follow from the fundamental 
fluctuation-dissipation theorem applied to EM interac- 
tions [22]. 

For isotropic fluctuations, the lowest order averaged 
change of the potential energy V{r) at location r is given 
by [211]: 



(A;^(r)), = i (Ar2),vV(r) 



(3) 



where (Ar^)s is the mean square deviation of the elec- 
tron at r executing a Brownian motion as driven by the 
random electric field satisfying. The resulting motion is 
divergent if the integration in Eq. (jlj is unregularized. 

However, the semiclassical coupling between the non- 
relativistic electron and the Fourier components of 



stochastic EM field becomes invalid as the corresponding 
frequency exceeds the Compton frequency of the parti- 
cle. Therefore the upper bound of integration is cut off 
near the Compton frequency luq = rnc^/H at aujc when 
coupled to a particle of mass m. Here a is an order one 
cutoff parameter that can be fixed either by experiment 
or renormalizable quantum calculation, yielding a « 3/4 
for the EM Lamb shift. This then leads to a finite overall 
energy shift through the expectation value: 



AE = ((Ay(r)),)* 



(4) 



where 
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{fir)h- /(r)|*(r)|^d' 



using the time- independent wavefunction ^'(r) of the 
quantum state. 

Position oscillations under spacetime fluctuations. — 
Despite the heuristicity of the semiclassical approach to 
the position oscillation of a non-relativistic electron in 
vacuum summarized above, it is supported by experi- 
ment and is justified in the subsequent full-fledged QED. 
To proceed, we wish to construct a gravitational analog 
that may one day be fully justified by a consistent, or 
sufficiently effective, theory of quantum gravity. To go a 
step further, we would like to outline a new experiment 
that could verify the resulting prediction and to hopefully 
guide theoretical development of quantum gravity. 

As a relatively simple, yet generic, model for spacetime 
fluctuations let us consider the perturbed metric g^i, = 
Vfj-f + hfj-v in the transverse traceless gauge of linearized 
gravity, where ii,v = 0, 1, 2, 3 and rj^^ = diag(— 1, 1,1,1) 
is the Minkowski metric and the metric perturbation h^i, 
is purely spatial, i.e. ho^ = 0. For i,j = 1,2,3, the 
metric perturbation hij may be interpreted as the gravi- 
tational radiation fleld. Up on applying the fundamental 
fluctuation-dissipation theorem in the fashion of , the 
thermal fluctuations of this fleld satisfies 
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E{uj,T) dco 



(5) 



where E{lu,T) is the same as given by Eq. (See e.g. 
23.) 



The effect of the fluctuating fleld hij is to induce a 
homogeneous and isotropic strain fluctuation of geom- 
etry that in turn induces a proper distance fluctuation 
between two free fall trajectories. If the two trajectories 
have a mean relative spatial displacement vector r with 
a constant norm r, then their relative position oscillation 
due to the fluctuating field hij can be described as an 
isotropic inhomogeneous fluctuation with a mean square 
deviation given by 



(Ar2)s 
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(6) 
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See Fig. [T]for illustration. In stark contrast to the ho- 
mogenous position oscillation induced by fluctuating EM 
field jU, the position oscillation induced by fiuctuating 
gravitational field increases with the mean separation. 

This is analogous to the increased effect of normal grav- 
itational waves on a detector using a pair of free masses 
with a large separation. For zero temperature T — 0, 
substituting Eq. ([5]) into Eq. ([6]), and imposing auJc as 
the high frequency cutoff following a similar argument as 
the EM case to ensure the consistency of the semiclassical 
treatment, we get 



Ar^ 3 = -i; jr^ 

671 TOp 



(7) 



where mp = y/hc/G is the Planck mass. Here the free 
cutoff parameter cr is also expected to be order one. 

Is there any gravitational Lamb shift in a free hydrogen 
atom? — Having established the mean square derivation 
for the position oscillation due to zero point fluctuations 
of spacetime (O , we are now in a position to apply it to 
physically interesting cases. A natural first case to study 
is to compare the gravitational Lame shift with the EM 
Lamb shift of the electron in a hydrogen atom using the 
same formula Q, but with different position oscillation 
characteristics. 

The potential energy in question takes the Coulomb 
form V{r) = —e^jr where r is the position of the electron 
relative to the CM of the atom with norm r, and e is the 
electron charge. Since V^l^(r) = 47re^(5^(r) it follows 
immediately from Eq. (jH) and Eq. Q that l\E = 0. 

Therefore, to the lowest order there is no gravitational 
Lamb shift for electron energy levels in a simple atom. 
The intuitive reason is that the position oscillation sim- 
ply vanishes at r = in this case, but V^F(r) vanishes 
everywhere except r = 0. It is instructive to compare 
V^F(r) to the material distribution of a gravitational 
wave detector. Analogously, one cannot detect gravita- 
tional waves with a free point mass alone! 

Gravitational Lamb shift of BEG in a compact trap. — 
To obtain a nonvarnishing gravitational Lamb shift, it is 
necessary to consider a non- Coulomb confining potential 
V{r) with V^y(r) ^ for r > 0, assuming r = is the 
classical center of mass of the system. 

A physically interesting case fulfilling this condition 
is the BEC in a harmonic trap. In this case, the CM 
motion of the condensate decouples from its internal de- 
grees of freedom [2j|. This CM wavefunction satisfies a 
single particle Schrodinger equation with mass m as the 
total mass of the condensate. The crucial property be- 
ing invoked here is that the total Hamiltonian i?tot of 
the trapped BEC is completely separable into two non- 
interacting part: iJtot — Hqm + -ffint, where Hcm is the 
Hamiltonian for the CM motion and Hint is that for in- 
ternal motions. To first approximation the total energy 
shift is the sum of the shift of the two parts. However, 



the mass term associated with internal modes is of the 
order of the individual atomic mass implying negligible 
energy shifts according to Eq. (O. 

Let us consider a 3-dimensional spherically symmetric 
harmonic potential V{r) = ^ muj'^r'^ with an isotropic 
trap frequency u}/2tt. Clearly we have a positive constant 
V^y(r) = "imuP which promises a positive gravitational 
Lamb shift. It then follows from Eq. ([7|) that the local 
shift of potential energy ([3]) yields 



Al/(r) 



37r TTlp 
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(8) 



where r — + y'^ + z^. Evidently, the higher shifts 
occurs further away from the center of the potential. 

It is straightforward to see from Eq. Q that, for 
the CM quantum state 5'(r) — (j") with the 

conventional 3D harmonic oscillator quantum numbers 
rij;, UyjUz = 0, 1, 2, • • • in the Cartesian basis, the gravi- 
tational Lamb shift is given by 



l^h.. 



(9) 



From this expression we see that the resulting gravita- 
tional Lamb shift scales quadratically with the total mass 
TO of the trapped BEC and increases with larger quan- 
tum numbers nx,ny,nz. The reason for the mass scal- 
ing is twofold: the potential energy increases with mass 
according to Eq. ([5]) for the same trap frequency, and 
cutoff for integrating the fluctuating gravitational fleld 
using Eq. ([5]) has also increased with a larger Compton 
frequency loq . The energy shift increases with the larger 
quantum number due to a more off-center expectation 
value for the CM position of the BEC, allowing higher 
levels of the fluctuating tidal force relative to the center 
of the trap. 

For a realistic ground state BEC in a harmonic trap 
with = 0, however, the gravitational Lamb 

shift is fantastically small, as the total trapped mass to < 
10^ GeV/c^ falls far short of the Planck mass top ^ 10^^ 
GeV/c^. This situation can be compared to detecting 
gravitational waves with a small continuum, where the 
effect could be calculated, but might not be measured. 

Gravitational Lamb shift of BEG in an extended 
trap. — In order to obtain a measurable effect, it is nec- 
essary to envisage an experiment with a longer baseline 
which allows the interactions with larger tidal forces as- 
sociated with the zero point gravitational waves. Con- 
sider now a cylindrically symmetric harmonic potential 
V{r) = i mw^(x^ -I- y^) with a transverse trap frequency 
w/27r in the (x, y)-plane and a zero longitudinal trap fre- 
quency along the z-axis. However, the atoms are "boxed" 
to move only freely between —L/2<z<L/2. In other 
words, the wavefunction is subject to zero boundary con- 
ditions at z = zLL/2. Thus we have W'^V{r) = 2muj'^ 
within the trap. 
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It then follows from Eq. ([7]) that the local shift of 
potential energy ([3]) yields 



(10) 



The idea is to allow the wavefunction to extend to 
macroscopically large value of along the z-axis, where 
the gravitational fluctuations become much larger rel- 
ative to the vicinity of center of the potential. The 
CM wavefunction in this case takes the form ^'(r) — 
(a;, (z) where = 0, 1, 2, • • • are the 

quantum numbers for the effective 2D harmonic oscil- 
lator wavefunction ^„^_„j^ (x, y) in the (x, ?/)-plane, and 
= 0, 1, 2, • • • is the quantum number for the effective 
boxed particle state {z) in the z-direction. 

The amount of shift is also calculated using dU to be 

l\E = X 

Dtt mp 



1^ ^ 
^ " mcj 12 



1 %^ 



.(11) 



For a trap with a long length-to-width aspect ratio, 
L 3> ^Jh/mjjj, the energy shift is dominated by the sec- 
ond term in PT|) . Furthermore, for T > uK, a BEC 
containing atoms satisfying 



En, — 



2mL2 



implies a large value Uz > 10 x VN x (L/m) ^ 1 for a 
sufficiently large value of L. It is however important to 
note that this restriction is not severe, as even at absolute 
zero temperature with = 0, it only affects the energy 
shift by an order one factor. In these limits, the energy 
shift reduces to 



AE 



2t2 



27tt ml 



(12) 



For example, assuming cr 1, if = 27r kHz and N 
is the number of trapped rubidium atoms, then the ratio 
of the energy shift to the difference of transverse energy 
levels is 



AE 
tvjj 



4 X l{r^^N^{L/iaf 



(13) 



For N = 10^ and L = 1 cm, the gravitational Lamb 
shift is vanishingly small, with AE/hw 4x 10~^°. With 
a relatively moderate improvement of atom number to 
N — 10® while keeping L = 1 cm, the gravitational Lamb 
shift yields AE/hio ^ 0.04%. For N = 10® and a further 
increase of the trap length to L = 5 cm, the gravita- 
tional Lamb shift has a potentially measurable value of 
AE/hoj « 1%. 

The discussion presented here is readily extended to 
other trap geometries. We expect that such traps will 



be available in the near future and the proposed experi- 
ment could potentially lead to the first observable effect 
of low energy quantum gravity. Furthermore the deter- 
mination of the free cutoff parameter a through the en- 
ergy shift experiment provides a stringent criterion for 
any viable theory of quantum gravity to reproduce its 
measured value. 
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